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State selective preparation and manipulation of discrete-level quantum systems such as atoms, 
molecules or quantum dots is the ultimate tool for many diverse fields such as laser control of chemical 
reactions, atom optics, high-precision metrology and quantum computing. Rabi oscillations are one 
of the simplest, yet potentially quite useful mechanisms for achieving such manipulation. Rabi theory 
establishes that in the two-level systems resonant drive leads to the periodic and complete population 
oscillations between the two system levels. In this paper an analytic optimization algorithm for 
producing Rabi-like oscillations in the general discrete many-level quantum system is presented. 



PACS numbers: 3.65.Sq 
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I. INTRODUCTION 



During the past 20 years a number of methods has been devised for state selective preparation and manipulation of 
discrete- level quantum systems However, simple population oscillations, induced by a resonant driving 

pulse have received negligible attention as a prospective population manipulation method. This might be attributed 
to two reasons. The first is that Rabi theory is based on the rotating wave approximation (RWA), and all attempts 
to generalize it without RWA (e.g. @,Qil3) are mathematically very involved. The second is that no attempt has 
been made to analytically generalize the original Rabi theory beyond two-level systems. 

In this paper an analytic extension of Rabi theory to transitions in many-level systems is presented. The aim is to 
'design' a driving pulse of the form: 

F{t) = Fq m{t) cos {u}{t) t) (1) 

by establishing analytical optimization relations between its parameters: maximum pulse amplitude , pulse envelope 
shape m{t), and time dependent carrier frequency to(t). The goal of this enterprize is twofold. The first is to achieve 
as complete as possible transfer of population between two selected states of the system. The second is to make 
this transfer as rapid as possible. These two requirements, however, are contradictory; population transfer can be 
accelerated by using a more intense drive, but at the same time a stronger drive increases involvement of remaining 
system levels in population dynamics and hence deteriorates population transfer between a selected pair of levels. 

In this paper it is shown how, for a pulse of arbitrary shape and duration m{t), the drive frequency can be 
analytically optimized to maximize the population transfer amplitude between selected two levels. In other words, 
Rabi oscillation theory is re-formulated for the case of a many-level system driven by an arbitrary modulated pulse. 



II. THEORETICAL ANALYSIS 



All calculations in this section are done in a system of units in which h = 1. 



A. Calculation setup 

A quantum system with N discrete stationary levels with energies Ei (i = 1, ...jN) is considered. The system is 
driven by a time dependent perturbation given in Eq. In the interaction picture, the dynamics of the system 

obeys the Schroedinger equation: 

|a(0 ^ -iV{tMt), (2) 

where a(t) is a vector of time-dependent expansion coefficients ai(t), aAr(t). The NxN Matrix V(i) describes 
interaction between the system and perturbation. Explicitly, its elements are given by: 

V,j{t) = :^^m(i)(e*'*'^("(*)-'^'^)* + e-^^-^('^(*)+'^'^) *). (3) 

fiij is transition moment between the i-th and the j-th levels induced by the perturbation. Sij — sign{Ei — Ej) and 
LOij = \Ei — Ej\ are respectively the sign and the magnitude of the resonant frequency for the transition between the 
i-th and the j-th level. 

The aim is to induce population transfer between two arbitrarily selected levels, designated by a and f3, directly 
coupled by the perturbation (i.e. such that iiap ^ 0). To simplify equations, the time variable t is re-scaled to r, with 
transformation between the two given by: 

dr = ?^m{t)dt. (4) 
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Then with following substitutions: 



Eq. ^ transforms into: 



where: 



x{t) 



-{Lu{t) -UJij) 



2 

Sy- H-Wij) 



2 

Ma/3 



t{T) 



dT 



a(T) = -iW(r)a(T), 



(5) 

(6) 

(7) 
(8) 

(9) 
(10) 



Initial conditions for the problem of selective population transfer comprise complete population initially (at t = 
T = 0) contained in only one of the selected levels, either a or (3. The other selected level, as well as all the remaining 
N-2 'perturbing' levels of the system are unpopulated at this time. 

Population evolution Iii{t) of the i-th level is determined from Iii{t) — |ai(i)p. 



B. Population oscillations in the two- level system - recapitulation 

Having in mind that — —sji, the explicit form of the general dynamical equation (O in a two-level system is: 

Under certain conditions (for a thorough discussion see eg. 0) this equation may be simplified by introducing the 
rotating wave approximation (RWA). Within the RWA, dynamical impact of complex exponentials e^^"''*^'^-'^^'^' is 
neglected and these may be eliminated from the equation. Hence, reduces to: 



d 






d7 


a/3(T)_ 





d 

d^ 



gifc,l3{T)x{T)' 

g-ifa,fiir)x{T) Q 



0/3 (t) 



(12) 



Finally, adjusting the perturbation frequency to the resonant value (i^^T) — ujafs) results in /q/3(t) = 0, and l|12|) 
reduces to: 



(13) 



As may be easily demonstrated by solving this simple equation, the resonant perturbation induces complete periodic 
transfer of the population between two levels. The time O (in units of r) or T (in units of t) required for a single 
population oscillation is determined from: 



d 




— —i 


1 






d^ 






1 




a/3(T) 



e 



dT = 



m{t)dt = TT 



(14) 



/O ^ ^0 

This is a well known result which forms the basis of 7r-pulse theory. 



C. Population oscillations in the three-level system 



As will be shown in subsequent sections, the whole analytical approach to the maximization of the population oscil- 
lation amplitude in a general many-level quantum system may be reduced to discussion of a three-level system. Along 
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with two 'selected' levels a and /3, the system now discussed contains one additional 'perturbing' level, designated 
with index p. The only requirements on the system internal structure are that iJ,af3, IJ^/Sp 7^ and /i^p — 0. While the 
first two requirements are necessary, the last one does not reduce the generality of the final results to any significant 
extent and it is introduced for calculational convenience. 
The three-level version of the dynamical equation ^ is: 



d_ 
d7 



ap(r) 








(15) 













_ap{T)_ 



The adiabatic approach and perturbing level population 



In order to solve equation p5|l two assumptions shall be made. The legitimacy of both of them may be checked 
retrospectively from the final solution. The first one is that the RWA may be applied for transition a ^ /3 so 
that exponentials containing 30/3 (t) may be dropped from equation (|15() . A simple way to verify this assumption's 
legitimacy may be done using Rabi profile plots, as described in j9|. The second is that the dynamical time-scale 
of solutions for aa^/sir) is much longer than that for ap{T). This enables one to regard a^(r) as a slowly changing 
parameter in the dynamical equation for ap{T). The solution to this equation hence may be obtained in terms of a 
parameter whose value needs not be known beforehand. 

Introducing the first of the assumptions into 115|) and reformulating it slightly, a set of two coupled differential 
equations is obtained: 



d 




= —i 


Qifcfi(r)x{Ty 






d7 






g-ifc0{r)x{T) Q 







aT 

Now consider just equation 117(1 . The formal solution to this equation is: 

ap(r) = -iRpp f {e'^P"'^^'^'''^'''^ + e~'3(^pi^'Mr'))ap{T')dT' 



ap{T), 



(16) 
(17) 

(18) 



This integral cannot be precisely evaluated until the exact form of the solution api^r) and optimized perturbation 
frequency w{t) - through which fppir) and gi3p{T) are defined (©, ©) - are known. Certainly, these are not 
known before the final solution of the whole optimization procedure is obtained. However, introduction of the second 
assumption enables one to evaluate the partial contribution to the integral in ((18|l from some interval tq < r' < t 
within which changes in all these functions are so insignificant that functions themselves may be approximated by a 
constant value. 

The aim of the optimization procedure is effectively to eliminate the dynamical impact of the perturbing level on 
population transfer between the two selected levels. If this is achieved, then dynamics of subsystem (a, 0) will be very 
similar to the dynamics of the pure two level system. Hence the population oscillation period will be about O = tt, 
which is then the dynamical time-scale for aa^p{T). As optimizing variations of the driving frequency w(t) are caused 
exclusively by the changes in perturbation envelope amplitude m{t{T)) it is transparent that dynamical time-scales for 
//3p(r) and gppir) are of the same order as that for m{t{T)). The dynamical time-scale of m{t{T)) must be of the same 
order as that of act.^(r) or longer, for otherwise not even a single complete population transfer would be achieved. 
Hence, in the interval tq < t' < t such that r — tq << tt, the optimized functions a/jlr'), fpp{T') and gi3p{T') may be 
considered constant, ajsir') w a/3(To), ffSpir') ~ fppiT) and (?/3p(t') w gf3p{To). Finally, in this interval x{t') may be 
approximated by x(t') « x{to) + m(f(ro)) ^ ^o)' several simple steps the following result is obtained: 



ap(T) « c(ro) - spp app 



m{t{T)) 

1 - App{T) 



1 



Sppj 



1 - A/3p(t) 

- SppAppir) 



-2iSf)j,Uj(T)t{T)\ ifi3piT)x{T) 



apir) 



(19) 
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OJaP + 




coir)- 


- ^aP 







where: 

^Pp " oTTT" — TTT ' ^^^> 

(21) 

A/3p(r) = ■ (22) 

The constant term c(to) includes both the ap(ro) and the integration constant obtained by inserting the lower limit 
value To into the solution. Its exact value cannot be determined since integration cannot be analytically stretched 
over the whole interval from r = to r = O. However, since ap(r = 0) = this constant is of the order of average 
value of remaining rapid complexly rotating expression, which is very nearly equal to zero. 

If detuning from resonance is assumed small, the perturbation frequency may be approximated by uj{t) w ujap 
so that A^p(t) = 0. The obtained solution (|19() then immediately provides an upper limit for the magnitude of 
perturbing level population: 

\Up{r)\<alp{l + \60p\f. (23) 
In general S^p << 1, so the whole bracket can be reduced to 1, and the maximum amplitude of the perturbing level's 

and the perturbation may be considered weak; if a^^ ~ 1 , the perturbation is very strong. This result may also be 
cast into a convenient quantitative form: to keep the 'leakage' of the population to the perturbing level p below a 
certain limiting value Mp, the greatest drive intensity which may be employed is roughly: 



population is roughly cr^^. Hence, this quantity may be regarded as a parameter determining the effective strength 
of the perturbation applied to a particular transition: if a'i << 1, then the dynamical impact of level p is negligible 



= | ^^^"^~^^^^ lVMp. (24) 

M/3p 



2. Optimized driving frequency 

When solution 1)19(1 is plugged into Eq. 1)16(1 a single closed dynamical equation is obtained for a two-level sub-system 
{a, (3): 
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~ —i 




0/3 (t) 





gifcl3iT)x{Ty 



with: 



1 — 5[3p/S.f3p{T) 

Now a transformation of the (a, (3) sub-system vector is sought: 



"6a (r)' 




aa{T) 


bp{T)_ 




ap(T)_ 



(25) 



(26) 



(27) 



(28) 



(29) 

As this equation is identical to ((13(1 . the corresponding solution would represent complete population transfer os- 
cillations between levels a and (3. Introducing transformation 1(27(1 into 1(25(1 . the following equation is obtained: 



with: 



A(r) 



Pi{r) 
P2{t) 



such that the transformed sub-system vector satisfies: 
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Mr) 




bp{T)_ 



(30) 
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If this is to be fitted to form (|29|l . the following conditions must be fulfilled: 



d 



= 0, 



which can be compactly written as: 



X(r) - ^p,{r) = 
fcp{T)x{T) - (pi(r) - P2{t)) = 



-(/„^(rMr)) = -x(r) 



0, 
0, 



(31) 

(32) 
(33) 



(34) 



Integration of the last equation yields the formal solution for A/3p(t) from which optimized perturbation frequency 
lo{t) may be extracted: 



A/3p(t) = {si3aSf}p)alp- 
+ (s/3aS/3p)CT^p 



1 — Spp 



x{t) Jo (1 - A^p(r'))(l - SppAppir')) 



m{t{T))dT' 



(35) 



x{t) 



„2iSfij,U}(T'')t(T') 

1 - A/3p(T') 



g-2is^pCi;(T')t(r') \ 



1 - (5/3pA;3p(T') 



Since A^p(t) is generally a small quantity, this equation may be solved iteratively, using A^p(t') = as the initial 
value. The contribution to the total A^p(r) from the second integral, containing rapidly rotating complex exponentials, 
may be shown to be minor compared to the one from the first, real integral. However, the very fact that there is 
an imaginary contribution to the optimized perturbation frequency indicates that the optimization procedure simply 
cannot completely annihilate the dynamical impact of the perturbing level. Nevertheless, as will be demonstrated in 
the following section, it may be done to a very good approximation. 

Transforming l|35() back to the original time coordinate t and keeping only the first integral yields the approximate 
recurrent solution for A^p(t): 



A/3p(t) = (s/3aS/3p)crfl (1 - 6f3p)- 



tJo {I - A0p{t')){l - Sf3pApp{t')) 



dt' 



(36) 



Introducing the zeroth-order approximation A^p(r') = into this equation, the analytic expression for the first-order 
optimized frequency is obtained: 



Uj{t) = LOaP + (W/3p - 0Jal3){si}aSpp)app{l - 5pp)- / {m{t')) dt' 



(37) 



For all but the strongest perturbations (i.e. such that cr^p ~ 1 or greater) higher-order corrections are not needed. 

Note that the frequency shift in (|36|l and H37(l may be either away from the perturbing line or towards it. Which 
case it will be depends on the relation between the energies of the three system levels: if level (3 has either the 
highest or the lowest total energy, so that both transitions [3 ^ a and (3 ^ p are energy-wise either 'downwards' 
(s/3a = spp = +1) or 'upwards' {spa = spp = —1) then the shift will be towards the perturbing line; on the other 
hand, if the energy of level f3 is in between the other two energies so that these two transitions are in the opposite 
directions {spa = —spp), the shift will be away from the perturbing line. 



D. Population oscillations in the many-level system 



The approach presented in Sec. II. 3 can be easily extended to include multiple perturbation levels. The perturbation 
now couples each of the selected levels to a certain number of perturbing levels, and each of perturbing levels' dynamics 
is calculated independently from all the others. This may be done as long as the perturbation due to each single 
perturbing level is kept reasonably small (gauged by standards of the three- level system). Accordingly, the many-level 
analog of Eq. (|25|l is: 
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(38) 
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FIG. 1: Pulse envelope profiles. Time t is on the abscissa and m(i) is on the ordinate. The dash-dotted line corresponds to a 
square pulse, the dashed line to a sine pulse and the dotted line to a sine squared pulse. The pulse is switched on at t = To 
and lasts until t = T. In all three cases, the maximum value of the perturbing field intensity field, Fq, achieved at time 
is the same. In the case of the three-level system, it is such that = 2.0 while in case of the many-level system crfot = 0.2. 



with: 

" 1 — S 

" 1 — s 

X0{t) = m{t{T))Y,S(3pCrfipRpp J- . . xw. . (40) 

^ (1 - A^p(rjj(l - df3pAi3p[T)) 

where all quantities are defined analogously to the ones in the previous section, and complex rotating terms have been 
eliminated. 

The many-level equivalent of (|34|l is: 

|:(/a/3(rMr)) =(x„(T)-xMr)), (41) 

and the first-order solution is: 

L0{t) = UJaP + Z]^l i^aq - (^ap){Saf3Saq)<jl^il ~ Saq) (42) 

+ Ep=i {i^Pp-^af)){spaSiip)alp{l-5pp)^j j^{rn{t')fdt'. 

The total perturbation intensity for the case of a many-level system may be estimated by considering the maximum 
total population of all perturbing levels: 

m n 

a^^t = Max{Iiq{t- To <t <T)) + J2 Max{Uq{t; Tq < t < T)) 

q=l p=l 

■m n 

9=1 P=l 

Hence, if af^^ << 1 the perturbation of the many- level system is small; otherwise it is large. In the following section 
it shall be demonstrated that H43|l provides not just a qualitative, but also an excellent quantitative criterion for 
determination of the impact of perturbing levels on population oscillation dynamics. 



III. NUMERICAL SIMULATIONS 



In this section numerical simulations of system dynamics for resonant (i.e. un-optimized) and optimized (deter- 
mined from (|36|l and (03)) perturbation frequencies are presented and compared. Several pulse envelope shapes are 
considered: square pulse (m(t) = 1), sine pulse (m(t) — sin{rtt)) and sine squared pulse ('m{t) — (sm(f2t))^) (Fig. 1). 
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FIG. 2: Optimized frequency plots for the three-level case. Time t is on the abscissa, the perturbation frequency uj{t) on the 
ordinate. The total pulse duration is 7.25 ns. Two straight solid lines indicate the two resonant frequencies of the system, oj^ja 
(upper) and ujpp (lower). The remaining three lines are optimized frequencies for the three types of pulse: the dash-dotted for 
a square pulse, the dashed line for a sine pulse and the dotted for a sine squared pulse. 



A. Three-level system 



First a simple three-level system is considered. In this case a full iterative solution for the optimizing driving 
frequency is easy to calculate from Eq. H35|) (with complex contributions neglected). System parameters have 
the following values {a.u. = atomic units): ujpa — 0.017671 a.u., SjSa = 1, M/3q = 0.073 a.u.; ujj3p = 0.017611 a.u., 
spp = —1, M/3p = 0.098 a.u.. These system parameters correspond to the three ro- vibrational levels of the HF molecule 
in the ground electronic state: a = {v = 0,j = 2,m — 0), f3 = {v = l,j = 1, m = 0), p = (w = 2,j = 2, m = 0). 
System parameters are such that the optimizing frequency shift is away from the perturbing line. In all cases, the 
total pulse duration T — Tq equals 7.25 ns. 

In order to present clearly the improvement that optimization of driving frequency induces in population transfer 
between the two selected levels, the perturbation strength in following examples is set to an extreme value: a'^^ = 2.0. 
Fig. 2 compares evolution of optimized frequency u}{t) with two resonant frequencies of the system, ojpa and ujpp. 
In Fig. 3 resonant and optimized population dynamics are shown for each of envelope shapes. The increase in the 
amplitude of the population transfer between the selected two levels is obvious. 



B. Many-level system 



As an example of a many-level system, the set of ro- vibrational states of the HF molecule in the electronic ground 
state is considered. The numerical model used for the calculation of the system dynamics includes 310 levels (31 
rotational x 10 vibrational). It is based on the HF internuclear potential data and electric dipole moment data from 
[To| and ^3 respectively. The targeted transition is (u = 1, j = 1, m = 0) — > (f = 0, j = 0, m = 0)). In all cases, the 
total pulse duration T — To equals 4.84 ns. 

It was demonstrated in the previous subsection that optimization indeed leads to improvement of the population 
transfer dynamics, even when the perturbation is very large. However, as in such conditions complete population 
transfer is unattainable, these examples were more of a qualitative nature from the standpoint of population transfer 
control. The many-level system considered now is more realistic than the previous three-level one and the focus is 
shifted to quantitative predictions. Hence, the employed drive intensity will be much smaller so that results can be 
directly applied to population transfer control. Pulse envelope shapes are the same as in the three-level case (see 
Fig. 1). Maximum amplitudes of electric field are likewise equal in all three cases, but now they are chosen so that 
a^g^ — 0.2. Since the perturbation is relatively small, the optimized frequency may be determined from the first-order 
approximate solution (|42|l . 

In Fig. 4 the optimized dynamics of two target system levels is shown for each of three pulse envelopes. In all 
cases two things should be noted. First, the general shape of the optimized dynamics of each of the two selected 
levels is fairly close to pure sinusoidal oscillations. This is more so, the smaller the perturbation strength parameter 
cr^o( is. However, the complete population transfer is again not achieved because a certain share of the population 
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FIG. 3: Comparison of the resonant and optimized population dynamics. Graphs on the left side present numerical solution 
to system dynamics for each of pulse types with a resonant perturbation applied, — uj/sa- Graphs on the right hand-side 
present optimized dynamics. The top row corresponds to a square pulse, the middle to a sine pulse and the bottom row to a 
sine square pulse. For the sake of clarity, only a (solid line) and (3 (dashed line) populations are plotted while p population 
is omitted. Although the optimization clearly does not produce clean two-level dynamics, the increase in the amplitude of 
population oscillation is nevertheless evident. 

unavoidably ends up in perturbing levels. Second, the actual instantaneous loss of population transfer is close to (and 
actually smaller than) m{t)(T^g^, as indicated by dotted line in each of the plots. This shows that af^^ indeed is a good 
quantitative (and not just qualitative) indicator of relative drive strength. The argument which led to relation H24() 
may hence again be applied to determine the maximum intensity of the driving radiation to be employed if population 
losses to perturbing levels are to be smaller than some pre-defined amount. 

IV. CONCLUSION 

As was stated in the introduction, the aim of this paper is to explore and refine the use of Rabi oscillations as a tool 
in selective population manipulation of complex discrete-level quantum systems. The main aims of such manipulation 
are as great as possible population transfer and at the same time as short as possible population transfer time. From 
the simple two-level theory it is well known that the an increase in drive intensity yields a reduction in population 
oscillation period. However, the same theory can neither fully disclose all the limitations of this result that arise from 
the complexity of the internal structure of a many-level system, nor can it handle the unavoidable loss of population 
to the rest of the system. Results presented in this paper fill this gap: they enable one to determine the maximum 
possible drive intensity (and hence the lower limit of time) with which oscillations of pre-selected amplitude (say 
99%) may be achieved, and at the same time to minimize unavoidable losses of the population to non-targeted system 
levels. Finally, the method of Rabi spectra (see 9]) presents a simple, yet useful conceptual supplement to the analysis 
presented in this paper. 

In order to achieve the quickest possible population transfer between two pre-selected levels, driving pulse should 
be tailored so that it produces only a single half-oscillation of the population. However, during research for this paper 
it has been noted that for strong fields standard 7r-pulse theory (see equation (|14() and reference is also deficient 
when it comes to the complex many-level systems. Work is currently in progress on analytical extension of standard 
TT-pulse theory that would resolve this issue. 
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FIG. 4: Optimized population dynamics of the two targeted levels in the many-level system. Time t is on the abscissa, 
populations na,/3(t) are on the ordinate. The top graph corresponds to a square pulse, the middle to a sine pulse and the 
bottom one to a sine squared pulse. In all cases a^^i = 0.2 and the total duration of the pulse is 4.84 ns. The dotted line on 
each graph indicates the function 1 — m{t)a'l^^. 
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